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Section  I.  GENERAL 
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1.  This  monthly  publication  summarizes  the  activities  of  the  Systems 
Analysis  Directorate.  The  purpose  of  this  note  is  to  give  wider  and 
more  timely  distribution  on  subjects  of  concern  to  the  conmand. 

2.  The  most  significant  Memoranda  for  Record  (MFR's)  and  other  techni- 
cal information  will  ie  published  as  notes  or  reports  at  a later  date. 

3.  In  order  to  assure  accurate  distribution  of  this  publication,  addi- 
tion or  deletion  of  addresses  to/from  the  DISTRIBUTION  LIST  are  Invited 
and  should  be  forwarded  to  the  address  below. 

4.  Inquiries  applicable  to  specific  items  of  interest  may  be  forwarded 

to  Commander,  US  Army  Armament  Materiel  Readiness  Command,  ATTN:  DRSAR-SA, 

Rock  Island,  IL  61201  (AUTOVON  793-4483/4628). 


DESCRIPTION  OF  THE  PROCESS  OF 
LASER  SPOT  MOTION  DURING  DESIGNATION 
OF  MOVING  TARGETS 


I I 


DEPARTMENT  OF  THE  ARMY 

HEADQUARTERS.  UNITED  STATES  ARMY  ARMAMENT  COMMAND 
ROCK  ISLAND.  ILLINOIS  61201 

REPLY  TO 
ATTENTION  OF: 

DPSAR-SA  £ 1 JAN  1977 

SUBJECT:  Description  of  the  Process  of  Laser  Spot  Motion  Daring 

Designation  of  Moving  Targets 


Commander 
Picatinny  Arsenal 

ATTN:  DRCPM-CAWS  (Messrs.  E.  Manley  & E.  Zimpo) 

Dover,  NJ  07801 


: 


1.  References: 

a.  MFR,  AMSAR-SAM,  16  Jan  76,  subject:  Laser  Guidance  Measurements 

Program. 

b.  MFR,  DRSAR-SAM,  20  Jan  77,  subject  as  above. 

2.  It  is  recognized  that  there  is  widespread  interest  in  the  principal 
sources  of  error  in  the  guidance  accuracy  of  Copperhead  and  HELLFIRE. 
Further,  it  has  been  established  that  laser  spot  motion  (jitter  or 
dither)  is  often  the  primary  external  error  source.  Over  the  past  year 
attention  has  been  given  to  measuring  the  pheonomenological  components 
of  spot  motion  (Ref  a.). 

3.  With  a recognition  of  the  importance  of  this  subject,  DRSAR-SAM 
has  prepared  a memorandum  which  examines  the  spot  motion  phenomenon 
in  greater  detail.  This  attached  memorandum  (Ref  b.)  is  submitted 
herewith  for  your  information.  The  MFR  should  be  of  interest  to 
systems  analysts  and  to  test  personnel  involved  in  the  laser 
measurements  propram. 

4.  Briefly,  the  MFR  treats  a statistical  technique  for  estimating  the 
dynamic  components  of  spot  mo  :ion,  makes  comparisons  between  an  analytic 
(theoretical)  model  and  experiment  in  both  the  frequency-  and  time-domains, 
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DRSAR-SA 

SUBJECT:  Description  of  the  Process  of  Laser  Spot  Motion  During 

Designation  of  Moving  Targets 

and,  finally,  derives  the  equations  used  in  a digital  computer 
simulation  of  spot  motion. 

FOR  THE  COMMANDER: 

s?na> 


1 Incl  M.  RHIAN 

as  Acting  Director 

Systems  Analysis  Directorate 

CF: 

PTA,  DRCPM-CAWS  (Messrs.  B.  Earrett  & J.  Williams) 

AVSCOM,  DRASV-WR  (Mr.  D.  Dunlap) 

MICOM,  DRSMI-RGT  (Mr.  Lewis) 

MICOM,  DRSMI-RGN  (Mr.  Pastrick) 

RSA,  DRCPM-CAWS-FO  (COL  Nulk) 

WSMR,  STEWS-TE-PC  (Steve  Kadner) 

AMSAA,  DRXSY-AAM  (Pat  Hill) 

AMSAA,  DRXSY-GS  (Dave  Barnhardt) 

DRSAR-SAM  (Mr.  Schlenker) 

DRSAR-SACF 
DRSAR-SA  RF 

✓Drsar-sam 

DRSAR-S/H  RF 
SARRI-RLR  (Dr.  Atuoruso) 

DRXHE-RI  (Mark  Brauer  & Dick  Carter) 


DRSAR-SAM 


20  Jan  1977 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Description  of  the  Process  of  Laser  Spot  Motion  During  Designation 

of  Moving  Targets 


1.  A set  of  references  related  to  the  above  subject  is  found  at  Attachment  1. 

2.  Background. 

In  connection  with  the  support  given  the  Copperhead  program,  the  under- 
signed has  repeatedly  emphasized  the  importance  of  laser  spot  motion  as 
a primary  error  source  in  guided  projectile  performance.  References  [1] 
and  [2]  illustrate  this  emphasis.  As  a consequence  of  the  importance 
attached  to  spot  motion  and  otner  aspects  of  the  laser  signature,  con- 
siderable efforts  have  been  sp^nt  studying  tracking  error  for  various 
designators  and  trackers.  Some  analytic  models  of  the  tracking  error 
are  given  in  References  [3]  and  [4].  In  particular,  Reference  [4]  describes 
the  statistical  techniques  used  to  produce  an  analytic  model  descriptive 
of  the  tracking  error  of  the  Ground  Laser  Locator  Designator  (GLLD)  observed 
during  the  CLGP  OT-1  Tracking  Tests.  All  of  these  models  are  suitable  for 
a digital  simulation  of  the  stochastic  process  of  spot  motion,  which  can 
be  used  as  input  to  terminal  guidance  simulations. 

3.  An  accurate  characterization  of  laser  spot  motion,  whatever  its  source, 
is  also  important  because  (stochastic)  process  dynamics  are  also  pertinent 
to  the  problem  of  laser  spillover  and  its  consequence  for  guidance  accuracy. 

In  Reference  [5],  which  addresses  this  problem,  spot  motion  dynamics  were, 
simply,  assumed. 

4.  In  an  effort  to  standardize  laser  designator  testing  procedures, 

increased  attention  (Reference  [6])  has  been  given  to  various  phenomena 
which  give  rise  to  laser  spot  motion:  device  jitter,  atmospheric  propagation 

effects  (scintillation),  and  designator  tracking  error.  In  a field-test 
environment,  these  phenomena  all  contribute  and  cannot  be  isolated  by 
direct  measurement.  Consequently,  some  of  the  past  analyses,  which  have 
considered  "spot  motion"  and  "tracking  error"  as  interchangeable  terms, 

have  been  somewhat  indiscriminate.  Pragmatically  this  lack  of  attention 
to  the  components  of  spot  motion  may  not  have  been  of  great  significance 
in  that  the  joint  process  is  what  is  of  concern  in  guided  projectile 
accuracy  analysis.  However,  a closer  attention  to  the  components  of  spot 
motion  does  produce  clarity  in  understanding  the  disparate  phenomena. 

5 . Outline  of  Memorandum . 

This  memorandum  reexamines  some  of  the  spot  motion  data  from  the  CLGP 
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SUBJECT:  Description  of  the  Process  of  Laser  Spot  Motion  During  Designation 

of  Moving  Targets 

OT-1  Tracking  Tests  with  the  intention  of  isolating  by  statistical  means 
the  component  attributable  to  the  human  tracking  error  from  those  produced 
by  other  mechanisms.  For  notational  simplicity,  throughout  the  balance  of 
this  paper,  device  jitter  and  beam-steering  effects  due  to  the  atmosphere 
will  be  subsummed  by  the  term  "scintillation  component."  This  over- 
simplification is  considered  justifiable  because  of  the.  dominance  of 
atmospheric  effects  in  the  examples  treated.  In  decomposing  laser  spot 
motion  into  tracking  and  scintillation  components,  we  exploit  the  fact 
that  these  components  have  different  dynamics  which  are  nearly  separable 
in  the  frequency  domain. 


6.  Because  the  analysis,  comparison  with  experimental  data,  and  algorithmic 
description  of  spot  motion  components  is  lengthy,  these  mathematical 
developments  are  divided  into  five  separate  and  nearly  independent 
annexes : 


Annex  1.  Decomposition  of  Laser  Spot  Motion  into  Tracking  and  Scintillation 
Components. 

Annex  2.  The  Analytical  Covariance  Function  for  Laser  Spot  Motion  Developed 
From  the  Spectral  Density. 

Annex  3.  Digital  Computer  Implementation  for  the  First-Order  Component 
of  the  Spot  Motion  Error. 

Annex  4.  Digital  Computer  Implementation  for  the  Second-Order  Component 
of  the  Spot  Motion  Error. 


Annex  5.  Autospectra  of  the  Digital  Error  Processes. 


7.  Summary . 

In  Annex  1 it  is  argued  that  the  autospectrum  of  spot  motion  consists 
mainly  of  separable,  additive  components:  one  due  to  the  human  tracking 

process,  possessing  approximately  second-order  dynamics,  and  the  other 
due  to  processes  having  approximately  first-order  dynamics  with  a signif icantly 
higher  crossover  frequency  than  the  human  tracker  component.  Expressions 
for  the  gain  constant  in  each  component  of  the  autospectrum  are  developed 
in  terms  of  the  variance  of  each  component.  A numerical  example  of  the 
analytic  autospectrum  is  compared  with  an  experimental  estimate  derived 
from  the  azimuthal  error  of  Run  75A  of  the  CLGP  OT-1  Tracking  Tests. 

Parameter  values  of  the  analytic  model  were  selected  by  trial  and  error, 
under  constraints,  to  produce  an  acceptable  Chebychev  fit  to  the  experimental 
autospectrum  in  the  sense  of  minimizing  the  maximum  residual. 


8.  It  is  sometimes  useful  to  rake  comparisons  between  analytic  models  and 
experimental  estimates  in  the  time  domain  as  well  as  in  the  frequency  domain. 
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SUBJECT:  Description  of  the  Process  of  Laser  Spot  Motion  During  Designation 

of  Moving  Targets 

One  place  where  such  a comparison  is  useful  is  in  guided  flight  simulations 
in  which  a long  experimental  record  of  spot  motion  is  employed  directly 
as  input.  With  the  autocorrelation  (or  autocovariance)  function  one  can 
maximize  the  number  of  statistically  independent  replications  obtainable 
from  the  record.  Annex  2 develops  an  expression  for  the  autocovariance 
function  for  the  (joint)  spot  motion  process  and  extends  the  example  of 
Annex  1 by  comparing  the  autocorrelation  function  of  the  analytic  model 
with  an  experimental  estimate. 

9.  Annex  3 treats  the  digital  implementation  of  the  first-order  component 

of  spot  motion.  Implementation  is  achieved  by  low-pass  filtering  of  gaussian 
white  noise.  Expressions  are  developed  for  the  coefficients  of  the  digital 
filter  and  for  the  variance  of  the  input  noise  required  to  produce  a given 
output  process  variance. 

10.  Annex  4 parallels  the  developments  of  Annex  3 while  treating  the 
second-order  component  of  the  spot  motion  process.  This  development  is 
the  basis  of  the  principal  algorithm*  which  DRSAR-SA  and  others  have  used 
to  describe  spot  motion  during  the  past  four  years. 

11.  Annex  5 develops  the  autospectra  of  the  scintillation  (first-order) 
and  tracking  (second-order)  dynamic  components  as  implemented  digitally. 

The  autospectrum  of  each  component  is  separately  compared  '/ith  its  con- 
tinuous-time (analog)  counterpart.  Some  sources  of  distortion  in  the 
digital  implementation  and  in  the  process  of  making  statistical  estimates  of 
autospectra  are  identified,  and  are  quantified  using  the  first-order 
process  as  an  example. 


GEORGE  J.  5CHLENKER 
Operations  Research  Analyst 
Methodology  Division 
Systems  Analysis  Directorate 


*A  special  case  is  presented  in  subroutine  D1TH  of  Reference  [3]. 
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ANNEX  1 


DECOMPOSITION  OF  LASER  SPOT  MOTION 
INTO  TRACKING  AND  SCINTILLATION  COMPONENTS 


Let  the  variance  of  the  spot  motion  at  the  target  due  to  only  tracker 
2 2 

error  be  a in  (m  ) . Experienc  e with  the  TOW  and  GLLD  trackers  as  well 
as  airborne  man-operated  trackers  shows  that  tracking  error  exhibits 
approximately  second-order  dynamics. 

Let  the  tracking  process  be  second-order  Butterworth  so  that  the 
spectral  density  is  given  by 


H“  (V)  = 


1 + (v/vt) 


, 0 < V <“  1 


with  gain  constant  A and  natural  corner  frequency  vt- 
The  variance  or  total  power  of  the  process  is 


oo 

O2  = H2 

c J0  C 


(v)  dV 


Whence 


w 

, f 

1 Jo  1 + 


= A V TT  /~2  / 4 


A = 4 a 


JT  TT 


A = 2 / 2 a2 


(m2/hz) 


(1-1) 


Normally,  the  corner  frequency  V is  determined  by  the  dynamics  of 

the  human  motor  system  and  lies  below  one  hertz. 

Scintillation  data  suggest  that  the  majority  of  scintillation  power 

lies  above  one  hertz  at  frequencies  to  which  the  human  cannot  respond. 

Consequently,  tracking  error  noise  and  scintillation  noise  are  substantially 

uncorrelated.  Data  on  intensity  of  scintillation  above  10  hertz  suggest 

this  process  has  nearly  first-order  dynamics. 

We  suppose  that  a first-order  process,  uncorrelated  with  the  tracking 

2 

process,  with  variance  0 is  superposed  thereon.  The  spectrum  for  this 

s 

scintillation  process  is  given  by 


H"  (V)  = 


1 + (v/\>g)  ‘ 


with  corner  frequency  V 


Therefore, 


H (V)  dv 
s 


O = Bv 
s s 


I — 

J0  1 + f 


2 „ TT 

a = B v — 
s s 2 


Whence, 


B «=  — — 

TT  V 

S 


(1.6) 


(1.7) 


(1.8) 


(1.9) 


(m  /hz). 


(1.10) 


(1-2) 
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If  the  processes  are  superposed,  the  spectrum  of  the  joint  process 

is  just  the  sum  of  the  component  spectra.  Letting  the  joint  autospectrum 

2 * 
be  H <'v), 

H2(V)  = H2  (V)  + H2  (v)  (1.11) 

L S 

H2(v) t + — r • (1.12) 

i + (v/v  r i + (Wvcr 

C 5 

Example 

Based  on  daytime  tracking  tests  at  WSMR  during  the  CLGP  OT  1 Tests, 
the  following  parameter  values  were  selected  from  Run  75A  for  the  azimuthal 
error : 

V,  * 0.7  hz  , V =3  hz_ 
t s 

a = 0.230  m , as  = 0.165  m (at  3 km  range) 

A - 6.8038  10~2  , B = 5.7773  10~3  (m2/hz) . 

Numerical  results  are  shown  in  Table  1 and  in  Figure  1. 

TABLE  1.  ANALYTIC  ESTIMATES  OF  THE  TRACKING  AND  SCINTILLATION 
COMPONENTS  OF  AZIMUTHAL  LASER  SPOT  MOTION  IN  THE  AUTOSPECTRUM 


V 

H2(v) 

H>) 

H2(V) 

(m2/hz) 

0.3 

6.582 

io~2 

5.720 

io’3 

7.154 

io"2 

0.5 

5.398 

io-2 

5.621 

io-3 

5.961 

1G~2 

0.7 

3.402 

lo"2 

5.479 

io"3 

3.950 

io-2 

1.0 

1.317 

10-2 

5.200 

io-3 

1.837 

io"2 

1.5 

3.081 

io-3 

4.622 

io-3 

7.702 

io-3 

2.0 

1.006 

io"3 

4.000 

io-3 

5.005 

io-3 

4.0 

6.375 

io-5 

2.080 

io-3 

2.144 

io“3 

An  alternative  notation  for  the  spectral  density  (autospectrum)  is  ^xx(v) » 
for  azimuth,  and  T (v) , for  elevation 
yy  (1-3) 
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Interval/r  (v)  (0.67 


Frequency,  v (hz) 


Figure  1.  Comparison  of  an  Analytic  Model  of  the  Autospectrum 
of  Azimuthal  Laser  Spot  Motion  With  an  Experimental  Estimate 


ANNEX  2 


THE  ANALYTICAL  COVARIANCE  FUNCTION  FOR 
LASER  SPOT  MOTION  DEVELOPED 
FROM  THE  SPECTRAL  DENSITY 


The  spectral  density  (autospectrum)  of  laser  spot  motion  in,  say, 
the  x-direction  due  to  the  joint  effects  of  tracking  jitter  and  atmospheric 
scintillation  (beam  steering)  is  given  by 


r (v)  = 


l+(v/v X l+(v/v  )' 

t S 


where 


2^2  of 


with  variance  contributed  by  tracking  Ot  and  variance  contributed  by 

2 

scintillation  0 . By  definition, 

s 


2 . 2 
o+o  = 

3 J0 


r (v)  dv. 

XX 


The  natural  corner  frequency  of  the  second-order  tracking  process  is 

V and  that  of  the  first-order  scintillation  process  is  v . 
t s 

An  alternative  expression  for  (2.1)  in  terms  of  the  angular  frequency  u 


is  I”  (ou) 

XX 


I”  (oj)  = r (v(oj))dv 

XX  XX 


(2-1) 
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r (to)  = ^r  (^-) 
xx  2tt  xx  2 u 


A U) 4 (2tt)~1  B to2  (2ir)_1 

r;x(w)  - -TV* — + — s-  2 

OJ  + (0  (0  + (0 

t s 


Now,  the  spectral  density  is  by  definition  the  Fourier  cosine  trans- 
form of  the  autocovariance  function  Y (t) 

xx 


r (u) 

xx 


> (t)  cos  tot  dt 

XX 


In  general,  the  Fourier  transform  pairs  f(t)  and  g(to)  are  related  by: 


f(t)  cos  (tot)  dt 


f(t)  = 2n  g(<o)  cot  (tot)  dto  . 


From  (2.5), 

OO 

y (t)  = tr  r vw)  cos  (tot)  dto  . 
xx  Jo  xx 

The  following  general  results  will  be  used  to  develop  an  expression 
for  Yxx(t)  for  the  spot  motion  process: 


OO 

Jo  ^ 


2 , 2.-1  — 1 —a  t 

L + to  ) cos  tot  dto  = (2  7T  a^)  e 1 


[,4.  4,-1 

i (a,  + o)  ) cos  tot  dto  = 

Jo  z 

. „ -a  t/v^  _ a t 

(2tt)  3 2 6 sin  ^ — ) • 


(2-2) 


Using  the  expression  for  f (to)  in  (2.4), 


i -to  t 

T (oo)  cos  tot  du)  ■ B v (2tt) ~ e 8 

XX  s 


_1  -dit/^2  to  t 

+ A v (2ir)  e sin  (—  + -~) 

t '4/2 


With  A and  B as  defined  in  (2,1)  and  with  (2.7), 


-to. 


t/J2 


xx 


O '*'*■  _ tl)  t 

Y (t)  = e sin  (^-  + -^-)  + Og2  e 8 


-to  t 


Th'i  autocorrelation  function  is  defined  as 


P (t)  * y (t)/y  (0) 

xx  'xx  xx 


In  this  case, 


Y (0) 

XX 


2 . 2 
a + a 
t s 


An  Example 

Using  parameter  estimates  from  tracking  run  number  75A  of  the  CLGP 


OT  1 Tracking  Tests 


■ vt = 

0.7  hz 

Wt  = 

4.398  r/s 

V - 

3 hz 

U)  = 

18.850  r/s 

s 

s 

o = 

0.23  m 

(at 

3 km  range) 

t 

o = 

0.165  m 

s 


a2  /(a2  + a2  ) = 0.6586 
t t s 

a2  /(a2  + a 2 ) = 0.3414 

S t S 

pxx(£)  = 0.9314  exp  (-3. lit)  sin  (0.7854  + 3. lit) 

+ 0.3414  exp  (~18.85t). 

This  result  is  compared  v/ith  -the  experimental  estimates  in  Figure  2 


I 


t* 


1 


i 


(2.14) 


( 2-4) 
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ANNEX  3 


1 


DIGITAL  COMPUTER  IMPLEMENTATION  FOR  THE 
FIRST-ORDER  COMPONENT  OF  THE  SPOT  MOTION  ERROR 

Marginal  statistics  for  the  first-order  process  are  assumed  gaussian. 

However,  for  brevity  the  word  "gaussian"  is  often  omitted  in  describing 
the  process.  One  implements  a first-order,  continuous  (or,  analog) 
stochastic  process  with  a first-order  digital  filter  having  a discrete,  mean-zero 
gaussian  white  noise  input.  The  design  of  the  filter  is  based  upon  the 
assumption  that  the  continuous-time  process  is  being  time  sampled. 

Continuous-time  (Analog)  Processes. 

The  following  first-order  transfer  function  H'(s): 


H*(s) 


(B/2ir)1/2  w 


s + U) 


(3.1) 


having  a continuous,  white  noise  input  generates  the  autospectrum 

2 


[H*  (<*»)] 


2 (B/2tt)  o.s 

“ 2,2 
to  + (1) 

s 


(3.2) 


as  in  equation  (2. A)  of  Annex  2,  where  the  autospectrum  of  the  output  is 
just  the  squared  modulus: 


[H * (o>) ] 2 * H'(j<o)  H’(-jto)  (3.3) 

for  a white  noise  input  of  unit  variance. 

For  convenience  a normalized  version  of  (3.1)  is  employed  such  that 
gain  is  unity  for  s equal  to  zero: 

H(s)  » to  /(s  + to  ) . (3.4) 

s s 

Attention  to  the  variance  of  the  input  white  noise  is  deferred  until  later. 
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Filter  Design. 

A time  sampling  of  the  continuous  process  characterized  by  (3.4)  can 
be  used  as  the  basis  of  the  digital  implementation.  The  z-transform 
describing  the  digital  implementation  of  (3.4)  can  be  obtained  by  performing 
a bilinear  transformation  from  the  s-plane  into  the  z-plane.  (See  Walsh 

[3.1]  or  Stanley  [3.2]).  The  bilinear  z-transform  is  taken  by  substituting 

s ■ (z  - l)/(z  + 1)  (3.5) 


into  (3.4). 

To  minimize  distortion  of  the  spectrum  of  the  digital  process  relative 
to  the  equivalent  analog  process  we  use  a warped  analog  cutoff  frequency, 
a (instead  of  Wg) , where 

a = tan(iogT/2)  = tan(ir  VgT)  ,**  (3.6) 

with  sampling  interval  T and  desired  angular  cutoff  frequency  U)g  or 
natural  cutoff  frequency  Vg. 

* 

Using  the  warped  cutoff  frequency,  the  digital  transfer  function  is 


V*> 


a 


z - 1 
z + 1 


+ a 


(3.7) 


The  digital  transfer  function  is  the  ratio  of  the  z-transforms  of  output 
to  input  of  a digital  filter. 

[3.1]  Walsh,  P.  J.,  A Study  of  Digital  Filters,  AD710381,  Naval  Postgraduate 
School,  Monterey,  CA. , Dec  1969. 

[3.2]  Stanley,  W.  D.  Digital  Signal  Processing,  Reston  Pub.  Co.,  Inc., 

Reston,  VA. , c.  1975. 

**This  is  a normalized  version  of  the  warping  transformation  of  frequency 
required  by  the  bilinear  transformation.  Normalization  is  achieved  by 
division  by  2/T. 
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The  useful  (cannonical)  form  of  (3.7)  is,  after  manipulation, 


Hz(z) 


aQ(l  + z-1) 
1 + bz-1 


(3.8a) 


a - a/ (a  + 1) 
o 


b “ (a  - l)/(a  + 1) , 


(3.8b) 


Remembering  that  z is  the  backspace  operator,  implementation  of  the 
desired  digital  filter  follows  directly  from  (3.8).  Notationally , let 
the  noise  input  to  the  filter  at  tine  index  t be  nt  and  the  output  at 
t be  xt,  with  t integer. 

Then, 


+ bx  , = a n + i n„  , 
t t-1  o t o t-1 


(3.9) 


x = a n + a n , - bxk  , 
t o t o t-1  t-1 


•(3.10) 


Since  the  desired  output  spectrum  requires  a white  noise  input. 


E[ntnt_kJ  =0  , k / 0 


(3.11) 


°n  • k ‘ °- 


Variance  Ratio  of  Output  to  Input. 

In  selecting  the  value  of  noise  variance  of  the  input  to  yield  a 
particular  output  variance,  it  is  necessary  to  know  the  ratio  of  output 
to  input  variance  of  the  filter.  In  the  following  development  an  expression 
for  this  ratio  will  be  derived.  It  will  be  useful  to  employ  the  following 
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notation  for  the  covariance  function: 


Y (k)  = E[x  x ] . 
'xx  t t-k 


In  this  notation,  we  desire 


(3.12) 


Y (°) /°  or  Yvv(0)/Y  (0). 
xx  n xx  nn 


Taking  the  mathematical  expectation  of  the  product  of  nt  with  both 
sides  of  (3.10)  yields 


since 


E[xtnt]  = a Qan  , (t  integer). 


E[n  n ] = 0 
1 t t-1 


(3.13) 


E[nt*t-il  = 0 • 


Similarly,  taking  the  expectation  of  the  product  of  nt_^  with  both  sides 
of  (3.10)  yields 


E[Vt-l]  = ao(1  " b)  an 


(3.14) 


After  squaring  both  sides  of  (3.10)  and  taking  the  mathematical  expectation, 


Y (0)  = E[x2]  = E[a2  n2  + 2 a2  n n , 
'xx  t 1 o t ot  t-1 

,22  ,.22  , 

+ ao  nt_]_  + b x^]  . 


2 a bn  xt  , - 2 a bn^  , x , 
o t t-1  o t-1  t-1 


(3.15) 


With  (3.13)  and  (3.14), 


wo)/wo)  - 2 v(1  + b)- 


(3.16) 
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Alternatively,  with  (3.8b)  , 


Yxx(0)/Ynn(0)  = a/(a  + • 


This  is  the  desired  result. 


(3.17) 


Summary. 


In  summary,  to  simulate  a first-order,  mean-zero  process  xt  having  a desired 
cutoff  frequency  and  standard  deviation  ag,  one  can  employ  the  digital 
filter  given  by  (3.10)  with  a gaussian,  white  noise  input  nt  having  standard 
deviation  given  by 


O = [Ynn(0)/Yvv(0)]  °« 

n nn  xx  s 


0 = [(a  + l)/a]  ' a , 

n s 


(3.18) 


with  a given  by  (3.6)  and  filter  constants  given  by  (3.8b). 
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ANNEX  4 


DIGITAL  COMPUTER  IMPLEMENTATION  FOR  THE 
SECOND-ORDER  COMPONENT  OF  THE  SPOT  MOTION  ERROR 


Analysis  of  a variety  of  man-operated  trackers,  such  as  those  reported 
in  [3],  indicates  that  the  marginal  probability  distribution  function  for 
the  tracking  component  of  laser  spot  motion  is  adequately  described  as 
gausslan.  Additional  experience  with  tracking  records  of  this  type 
indicates  that  a second-order  dynamical  system  characterizes  the  tracking 
error.  For  most  human  trackers  a good  digital  simulation  of  the  mean-zero 
portion  of  the  tracking  error  is  obtained  by  filtering  gaussian  white 
noise  with  a second-order,  low-pass  Butterworth  filter.  As  in  Annex  3, 
the  design  of  the  digital  filter  is  based  upon  the  assumption  that  the 
continuous-time  (analog)  process  is  being  time  sampled. 

Continuous-time  (Analog)  Processes. 

The  following  transfer  function  describes  a second-otuer  analog  process 
with  analog  corner  frequency  u)  and  damping  constant  C: 

cl 

CO2 

H (s)  m — T -2  (4.1) 

s + 2£  u)  s + w • 
a a 

Specifically  for  low-pass  But.erworth  filters  t,  = 1/^2  and 


Ht(s) 


32  + Jl  U)  S + OJ2 


(4.2) 


The  digital  transfer  function  associated  with  (4.2)  is  created  by 
mapping  from  the  s-plane  into  the  z-plane  using  the  bilinear  z-transform. 


The  digital  transfer  function  is  the  ratio  of  the  z-transforms  of  output 
to  input  sequences  of  a digital  filter. 
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(See  Walsh,  P.  J.,  Op  Clt  [3.1]).  The  bilinear  z-transform  is  taken  by 
substituting 


s = (z  - l)/(z  + 1) 
in  (4.2).  Thus, 


■ Ht 


H (z)  = 


2 2 
<*£  (z  + ir 

(z  - l)2  + /2  w (z  - 1)  (z  + 1)  + CJ2  (z  + 1)‘ 
a a 


Reduction  of  H (z)  to  the  form 
z 


-1  , -2 
a + a, z + a„z 

0 1 2 

-1  -2 

1 + b^z  + b2z 


implies  the  digital  filter,  since  z is  equivalent  to  a unit  backspace 
operator.  The  coefficients  in  (4.5)  are  given  by 


a = a.  = a)  /D 
o 2 a 


al  “ 2 ao 


b,  = 2(d)  - 1)/D 

1 a 


b - (U)Z  - Jl  to  + 1)/D 
2 a a 


D = d)  + *^2  d)  + 1 . 

a a 


Notitionally,  let  the  noise  input  to  the  filter  at  time  index  t be  nt 
and  the  output  at  t be  xt,  with  t an  integer.  Then,  from  (4.5), 

Xt  + Vt-l  + Vt-2  = Vt  + al"t-l  + Vt-2  ‘ (4*7) 
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Since  the  analog  angular  frequency  to  associated  with  the  bilinear 

n 

transformation  is  "warped"  or  distorted  relative  to  the  desired  digital 
cutoff  angular  frequency  we  employ  the  relationship  between  these 
parameters:  1 

w = tan(a)  T/2) 
a t 

or 

U)a  = tan(TTVtT),  (4.8) 

where  T is  the  sampling  period. 

Parenthetically,  note  that  a change  in  the  sampling  period,  for  a 
fixed  digital  cutoff  frequency,  would  change  to  and,  by  (4.6),  would 

cl 

change  the  coefficients  of  the  digital  filter. 

In  operation,  the  filter  arithmetic  is  performed  as  follows: 


Vt + ai\-i + Vt-i  - Vt-i 


Tt-2 


with  t = 3,  4,  ... 

The  filter  is  initialized  by  assigning  = 0,  their  mean  values;  and, 

then,  cycling  through  a sufficient  set  of  inputs  to  remove  the  effect  of 
the  initialization  transient. 

The  form  of  (4.9)  is  referred  to  in  the  statistical  literature, 

* 

eg  [4.1]  and  [4.2]  as  a mixed  autoregressive,  moving  average  model  since 
the  output  xt  depends  upon  past  values  of  the  output  — xt  ^ and  xt  ^ — as 
* 

[4.1]  Box,  G.  E.  P.  and  Jenkins,  G.  M.,  Time  Series  Analysis:  Forecasting 

and  Control,  Holden-Day,  San  Francisco,  c.  1970. 

[4.2]  Jenkins,  G.  M.  and  Watt;,  D.  G. , Spectral  Analysis  And  Its  Applications, 
Holden-Day,  San  Francisco,  c.  1969. 
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well  as  upon  (moving  average)  terras  in  the  input  — nf  and 

For  the  required  autospectrum  of  xt , n£  must  be  white.  Thus, 

E[ntnt_k]  = 0,  k ^ 0 

= o^,  k = 0.  (4.10) 

Furthermore,  future  values  of  the  input  are  uncorrelated  with  the  present 
value  of  the  output,  ie, 

E[nt+kxt 1 = °»  k>0*  (4.11) 

Ratio  of  Output  to  Input  Variances. 

In  selecting  the  value  of  noise  (input)  variance  to  yield  a particular 
output  variance,  it  is  necessary  to  know  the  ratio  of  output  to  input 
variances  of  the  filter.  An  expression  for  this  ratio  will  be  derived. 

We  employ  a notation  used  by  Jenkins  [4.20]  for,  respectively,  the  auto- 
covariance and  crosscovariance  functions: 


Y,%<k)  " Elxtxt-k] 


Yxn(k)  * E|Vt-k'- 


In  this  notation,  one  desires 


(4.12) 


Y»<0)/%  Y**<0)/Ynn«»- 


This  result  is  obtained  from  (4.7)  by  successive  multiplication  of  both 
sides  by  xfc  nt,  and  nfc  and,  then,  by  taking  the  mathematical 
expectation  on  both  sides  in  all  equations.  This  produces: 


e; 


7; ! 


Y (1)  = (a  - b.a  )o*  . 

xn  1 x o n 


Thus, 


Vxx<°>  + (1  + b2>  Yxx(1> 


[a  a + a (a  - b a )]  O 
ol  Z 1 1 o n 


(A. 13) 


After  squaring  both  sides  of  (4.9)  and  taking  expected  values,  one 
obtains 


<!  - ' »2>  V°>  - 2 blb2  Yxx(1>  ‘ 


2 2 2 2 
a [-2a  a.b,  - 2a  a b_  - 2a»b, (a,  - b,a  ) + a + a,  + a_], 
n oil  o22  211  lo  o 1 Z 


(4.14) 


One  can  solve  (4.13)  and  (4.14)  for  the  desired  variance  ratio. 

C^B  - B C2 

\x(0)/W0)  = A.B-—  A,  * 


(4.15a) 


12  12 


with 


A1  = bl 


A2  " 1 - bj  - b’ 


B1  = ! + b2 


B2  " -2  blb2 


C1  = Vl  + a2(al  " hlao) 


C-  = -2a  a.b,  - 2a  a,b.  - 2a_b  (a.  - b.a  ) + a2  + a2  + a2.  (4.15b) 

2 oil  oz2  211  lo  o 1 Z 


For  numerical  stability  the  denominator  of  the  r.h.s.  of  (4.15a)  must, 
of  course,  be  non-zero.  This  implies  that  the  parameters  b^  and  b must 
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lie  within  the  region  defined  by: 


|b,|  <1 


bl  - b2  <X 


bi  + b2  > -1. 


This  requirement  is  equivalent  to  the  statement  that  the  poles  of  (4.5) 
lie  outside  the  unit  circle.  See  Box,  Op  Cit,  [4.1].  In  practice,  stability 
requirements  are  satisfied  for  the  problem  of  modeling  human  tracker  error 
using  sampling  rates  of  10  or  20  hertz. 

Summary. 

To  simulate  the  second-order,  mean-zero  component  of  a stochastic 
process  having  the  desired  cutoff  frequency  and  standard  deviation 
one  can  employ  the  digital  filter  given  by  (4.9)  with  a mean-zero,  gaussian 
white-noise  input  n^.  having  a standard  deviation  give"  by 


°n  " ^nn(0)/W°)]  °t  * 


(4.16) 


with  variance  ratio  given  by  (4.15)  and  with  filter  coefficients  given 
by  (4.6). 
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ANNEX  5 

AUTOSPECTRA  OF  THE  DIGITAL  ERROR  PROCESSES 


First  Order  Component 

From  the  digital  transfer  function  of  the  first-order  process  (equation 
(3.8a)), 


H (z) 
z 


ao(1+z"1) 


-1 


(5.1) 


1 + bz 

and  with  the  following  mapping  into  the  frequency  domain: 


z ■ e 
the  squared  modulus 


Ju)T  J2ttvT 


(5.2) 


H?(v)  - H (e22™1)  H (e'-l2™1) 
1 z z 


is 


_ 2 a (1  + cos  2irvT) 

H?(v)  = 7~ T2 


(5.3) 


1 + b + 2b  cos  2ttvT 
with  V defined  over  positive  frequencies  from  zero  to  the  Nyquist  frequency 


Vf=l/(2T). 


The  results  in  (5.1)  and  (5.3)  were  derived  from  a normalized  analog 


transfer  function  having  unity  gain  at  frequency  zero.  Consequently, 


H^(0)  ■ 1.  Then,  the  autospectrum  of  the  first-order  digital  component. 


Txx(v),  is  given  by 


rxx<V>  " HI(V>  rnn(v)  ’ 


(5.4) 


where  ^nn(\>)  Is  the  autospectrum  of  the  digital  white  noise  (input)  process, 


given  by 
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(5.5) 


r (V)  - v:1  , 0 < V < V , 

nn  in  i 


or,  from  (3.18), 


T (v)  - 2T  o (.i+l)/a 
nn  s 


r (v)  - LT(1+L>.  a2  o < v < v . < 

nn  » 2 s — — f 

c a 

o 

The  noise  spectrum  T is  comparable  to  the  gain  constant  B in  (1.10). 

nn 

As  T approaches  zero  approaches  B. 

Example 


It  is  interesting  to  compare  the  value  of  the  first-order  component 
autospectrum  for  the  analog  process  (equation  (1.6))  with  that  for  the 
equivalent  digital  implementation.  To  this  end,  take  the  parameters 
of  the  example  given  previously  in  Annex  1:  Vg_3  hz,  O8*0.165  m.  Assuming 

that  the  sampling  rate  (1/T)  of  the  analog  process  is  20  uz, 

T ■ 0.05  sec,  and,  by  (3.6), 
a - 0.509525. 


From  (3.8b), 


a - 0.337540 
o 

b - -0.324920  . 


Then,  (5.3)  becomes 


hJ(v) 


0.227866  (1  + cos  0.314159v) 

1.1055728  - 0.649839  cos  0 314159V 


(5.7) 


(5-2) 


This  result  is  compared  with  that  of  the  corresponding  analog  process 
in  Table  2.  Additionally,  squared  moduli!  for  other  sampling  rates  are 
displayed  in  Table  2. 


TABLE  2.  COMPARISON  OF  THE  ANALOG  AND  DIGITAL  SQUARED 
MODULI I OF  THE  FIRST-ORDER  ERROR  PROCESS 


V 

(hz) 

hJ(v) 

2 

Hx (v) , digital  at  sampling  rate: 

analog,  Br-1 

20  hz 

40  hz 

80  hz 

0.0 

1.0000 

1.0000 

1.0000 

1.0000 

0.3 

0.9901 

0.9915 

0.9904 

0.9902 

0.5 

0.9730 

0.9767 

0.9739 

0.9732 

1.0 

0.9000 

0.9119 

0.9030 

0.9007 

1.5 

0.8000 

0.8183 

0.8045 

0.8011 

2.0 

0.6923 

0.7109 

0.6968 

0.6934 

3.0 

0.500C 

0.5000 

0.5C0C 

0.5000 

4.0 

0.3600 

0.3297 

0.3531 

0.3583 

6.0 

0.2000 

0.1205 

0.1817 

0.1955 

8.0 

0.1233 

0.0267 

0.0984 

0.1171 

10.0 

0.0826 

0.0000 

0.0545 

0.0755 

r (V)  (m2/hz) 
nn 

5.777  10-3 

8.066  10-3 

7.031  10~3 

6.431  10"3 

(5-3) 
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As  noted  from  Table  2,  the  squared  modulii  of  ae  analog  and  digital 
processes  are  In  excellent  agreement  over  the  interval  (0  V < v^/2). 

For  frequencies  within  the  interval  (vf/2  < V <_  Vf)  and,  particularly 
at  the  upper  end  of  this  interval,  the  autospectrum  of  the  digital  process 
departs  significantly  from  that  of  the  corresponding  analog.  Since  there 
is  no  spectral  content  for  the  digital  realizations  above  v^,  ie,  all 
of  the  variance  of  the  digital  process  must  occur  at  frequencies  below 
v^,  the  values  of  the  digital  autospectrura  exceed  those  of  the  analog 
having  the  same  variance  for  low  frequencies.  Thus,  the  autospectrum 
of  a digital  process  produced  by  a low  sampling  rate  (and  low  v^)  will 
display  substantial  distortion  relative  to  that  of  the  corresponding 
analog  process  . This  point  is  illustrated  in  Figure  3 in  which  different 
digital  autospectra  of  the  first-order  stochastic  process  are  compared 
to  a corresponding  analog  au^ospectrum. 

Of  course,  if  the  contribution  of  the  first-order  process  to  the 
total  spot  motion  variance  is  small,  this  distortion  may  be  acceptable. 
However,  if  fidelity  to  the  corresponding  analog  process  is  desired, 
a general  rule  might  be  that  the  sampling  rate  used  in  generating  the 
process  digitally  should  be  20  to  30  times  the  coiner  frequency. 

Naturally,  the  dynamical  response  of  the  system  accepting  this  first-order 
noise  is  also  a consideration  in  selecting  a sampling  generation  rate. 


i 
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0 Digital,  sampling  generation  at  20  hz 

j ' T I I * i'll.  I > i 1 I il  l I.  ..I  I I 1 I 

^ Digital , sampling  generation  at  80  hz 
| Common  standard  deviation  0.165m 


Frequency,  v (hz) 

Comparison  of  the  First-Order  Autospectra  for  the  Analog  and  Corresponding 


Figure  3 
Digital  Processes 


I 


; 

| 

■ 

Generally  the  dynam  c components  of  a signal  such  as  laser  spot  position 
are  viewed  only  at  disc  ete  points  in  time,  namely  when  the  pulsed  laser 
flashes.  Thus,  any  spo  motion  due  to  a fundamentally  continuous  (analog) 
process  is  Inherently  t me  sampled.  Further,  this  sampling  interval 
T*  may  be  different  from  that  used  in  the  digital  implementation,  T.  To 
quantitatively  assess  what  effect  the  digital  generation  Interval  T has 
on  the  output  series  obtained  by  sampling  at  a time  Interval  T*,  one  must 
exploit  some  sampling  theory.  This  theory  is  developed  below  fot  an 
ideal  (delta-function)  sampler  and  applied  to  the  first-order  component 
of  the  stochastic  process. 


1 

! 


>•  : 
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The  Transfer  Function  ?f  a Time-Sampled  Process 

Whenever  a stochastic  process  is  sampled  (or  subsampled)  the  autospectrum 
of  the  sampled  process  may  be  substantially  altered  relative  to  the  original 
process.  This  distortion  will  occur  whenever  the  original  process  has 
an  appreciable  variance  (or  power)  invested  in  frequencies  (v)  above  the 
Nyquist  (or  folding)  frequency  of  the  sampler.  If  the  sampler  is  operating 
at  a sampling  rate  (T*)~l,  the  folding  frequency  will  be  Vf  = (2T*)“* 
and  variance  in  the  original  signal  associated  with  V > Vf  will  be  con- 
founded with  the  variance  associated  with  v for  0 <_  v < Vf . In  this 
section  a theorem  for  ideal  (point)  samplers  will  be  applied  to  several 
processes  and  examples  will  be  offered  to  quantify  the  distortion  accompanying 
sampling.  In  doing  this,  the  author  follows  the  notation  and  results  of 
Kuo  and  Kaiser  [5.1],  p.  222  ff. 

Notationally , let  f (nT*)  be  the  value  of  the  original  process  evaluated 
at  points  in  time:  t = nT*,  with  n an  integer  and  T*  the  constant  sampling 
interval.  Further,  let  F(s)  be  the  Laplace  transform  of  the  original 
process  and  let  F*(s)  be  the  Laplace  transform  of  the  sampled  process,  ie,  of 
the  discrete  series  result  ng  from  sampling.  Then,  for  zero  initial 
conditions, 

F*(s)  = E”=_roF( s H jn27r/T*)  , (5.8) 

with  s “ Jio  defined  for  < w < °°  . 


(5.9a) 


P 

y 

■ 

: 

And, 

F*(jto)  = F(jw)  + E”=1  F(J(wHon))  + Z~ml  F(J(u-un))  , 

with 

u «=  2irn(T*)_1  (5.9b) 

n 

and 

o>  = 2ttv  . (5.9c) 

Since  the  modulus  of  F(jaj)  will  generally  decline  with  co  for  values 
of  cu  greater  than  some  value,  say,  tu* , it  will  suffice  to  approximate  F* 
with  a finite  (and  reasonably  small)  number  of  terms  in  the  infinite  sums. 
Thus, 

F*(jw)  a F(jw)  + Z“*  F(j(orfw  ))  + En*  F(J(u-u  )).  (5.10) 

n3i  n n=i  n 

, In  evaluating  the  autospectrum  of  the  sampled  process  it  will  be 

convenient  to  separately  find  the  rtal  and  imaginary  parts  of  each  of  the 
terms  in  F*(joj). 

Then, 

Re  (F*(ju))}  s Re  (F(jw)} 

l 1 

+ In*  Re  {F(j  (a  Ho  ))} 
n=l  n 

! i 

+ En*  Re  {F(j  (( >— a)  ))}  , (5.11) 

n=i  n 


and  likewise  for  the  imaginary  part  of  F*(jtu),  Im  {F*(jcj)}. 

Finally,  the  autospectrum  of  the  sampled  process  is  proportional  to 


, 


- Ira  (F*(ju))}  = Woj  ) [(oj/co  )2  + l]"1 

s s 

+ E . ((w+ti)  )/o)  ) [(urt-uj  )2/co2  + 1]  1 
n=l  ns  n s 

+ E , ((cj-u)  )/o>  ) [ (ux-co  ) /to2  + 1]  1 . (5.16) 

n=l  ns  s s 

This  case  has  been  e^ aluated  for  the  following  specific  parameters: 

V ■ 1,  3 hz  and  T*  = 0.1,  0.05,  0.025  sec.  The  results  are  dispayed  in 
Figures  4 and  5. 

A second  example  of  the  effect  of  sampling  is  the  case  in  which  a 
digital  implementation  of  a first-order  process  is  subsampled,  ie,  in  which 
the  sampling  frequency  (T*)  ^ is  a submultiple  of  the  digital  generation 
frequency  T The  digital  process  with  transfer  function  ^(z)  in  equation 
(3.8a)  will  be  the  example  taken  here! 

a (1  + z *) 

H (z)  = — r . (5.17) 

2 1 + bz'1 

With  the  mapping  into  the  s-plane: 


m 


sT 

* c 


» 


z 


(5-9) 

47 


(5.18) 


analog 


T sampled 


Frequency,  v (hz) 

Figure  4.  Comparison  of  a First-Order  Analog  Squared  Modulus  With  the  Squared 
Modulus  of  the  Ideal  Time-Sampled  Series  (Parameters:  v « 1 hz,  T*  *•  0.1  sec) 
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Figure  5.  Comparison  of  a First-Order  Analog  Squared  Modulus  With  the  Squared  Modulus 
of  the  Ideal  Time-Sampled  Series 

(Parameters:  v - 3 hz,  T*  » 0.05,  0.1  sec)  ' ' 


F(s)  = H (e8T) 
z 


aQ(l  + e"sT) 


1 + b e 


-sT 


And, 


F(ja>) 


aQ(l  + exp(-jWT)) 
1 + b exp(-juT) 


with  (from  (3.8b)) 


aQ  * a/(a  + 1) 


(a  - 1)  / (a  + 1) 


a *=  tan  (tp^T) 


In  this  case. 


a (1  + b)(l  + cos  U)T) 
Re  { F ( jto)  } - ° 


1 + b + 2 b cos  toT 


and 


a (1  - b)  sin  uT 
-Im  (F(ju)}  - ° 


(5.19) 


(5.20®) 


(5.20b) 


(5.21a) 


(5.21b) 


1 + b + 2 b cos  o)T 
foe  - tt  £ ujT  <_  it  ; and  when  wT  is  outside  of  this  region,  F(jw)  is 
identically  zero. 

Obviously,  the  sampling  interval  T*  must  be  an  integral  multiple  of 
the  generation  interval  T: 

T*  = mT  , m = 1 , 2 , 3 . . . 

In  this  example  F*(jio)  is  evaluated  for  Vg  = 3 hz  and  T » 0.025  and 

T •=  0.05  sec  with  m = 2.  Results  are  plotted  in  Figures  6 and  7. 

Thu  autospectra  for  the  time-sampled  process  in  Figure  7 may  be  compared 
with  the  unsampled  autospectra  shown  in  Figure  3.  The  effect  of  sub- 
sampling order  (m)  on  the  squared  modulus  is  shown  in  Table  3. 
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analog,  w/o  sampling 


digital,  sampling 
rate  40  hz 
T = 0.025  sec 


digital,  sampling 
rate  20  hz 
T*  = 0.05  sec 


Frequency,  V (hz) 

Figure  6.  Comparison  of  Squared  Modulii  for  Subsampled  Digitally  Generated 
First-Order  Processes  Having  Different  Sampling  Rates 
(Parameters:  v = 3 hz,  T *=  T*/2)  ,, 
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Figure  7.  Comparison  of  the  Autospectra  of  Time-Sampled  Processes:  Analog  and 
Digital  Implementations  (Parameters:  v c 3 hz,  Og  «=  0.165  m,  T*  = 0.05  sec) 


TABLE  3.  EFFECT  OF  SUBSAMPLING  ON  THE  SQUARED 
MODULUS  OF  A DIGITAL  IMPLEMENTATION  OF  A 
FIRST-ORDER  DYNAMIC  PROCESS  SAMPLED  AT  20  HZ 

Parameters:  V = 3 hz,  T*  ■»  0.05  sec 


digital  |F(jaj)  | 2 with  T = T*/m 
subsampline  order  m 


0.9987 

0.9949 

0.9887 

0.9800 

0.9690 

0.9560 

0.9409 

0.9240 

0.8855 

0.8420 

0.7709 

0.6475 

0.5306 

0.4276 

0.3406 

0.2687 

0.2100 

0.1624 


2 

3 

0.9987 

1.0311 

0.9950 

1.0272 

0.9887 

1.0208 

0.9801 

1.0120 

0.9692 

1.0008 

0.9562 

0.9875 

0.9412 

0.9721 

0.9440 

0.9549 

0.8861 

0.9157 

0.8430 

0.8714 

0.7724 

0.7991 

0.6500 

0.6737 

0.5342 

0.5551 

0.4323 

0.4507 

0.3463 

0.3625 

0.2753 

0.2898 

0.2174 

0.2305 

0.1704 

0.1824 

1.0537 

1.0497 

1.0432 

1.0342 

1.0229 

1.0093 

0.9937 

0.9762 

0.9363 

0.8913 

0.8177 

0.6901 

0.5694 

0.4632 

0.3735 

0.2995 

0.2392 

0.1902 


I 


Aside  from  the  distortions  in  the  autospectrura  (and  process  dynamics) 
due  to  digital  sampling,  there  are  distortions  in  the  spectral  estimates 
created  by  the  statistical  techniques  employed  in  analyzing  time  series 
data.  These  distortions  are  discussed  below. 

The  Effect  of  Lag  Window  in  the  Estimated  Autospectrum 

In  estimating  an  autospectrum  from  time  series  data,  it  is  essential 
to  provide  a means  of  aviraging  or  smoothing  spectral  estimates  to  insure 
stochastic  convergence  ai  the  length  of  the  series  (or  record)  grows 
indefinably  large.  For  a single  record,  x(t),  the  smoothing  is  efficiently 
done  by  applying  a weight  function  or  "lag  window",  w(t),  to  the  estimated 

A 

autocovari3nce  function,  y (t),  before  taking  the  Fourier  transform  to 

A 

form  the  smoothed  autospectrum,  r (id) . See  Jenkins  and  Watts,  1968,  Op  Cit. 
Thus,  for  a continuous  record  of  length  t e» 


t 

P (id)  = — y (t)  w(t)  cos  cot  dt.  (5.22) 

xx  ir  I xx 

Without  applying  a weight  function  to  the  theoretical  autocovariance, 
and  for  an  essentially  infinite  record,  one  would  obtain  the  theoretical 
(unsmoothed)  autospectrum: 


r (id) 
xx 


oo 

ll 

" Jo 


Y (t)  cos  (dt  dt. 

XX 


(5.23) 


Parenthetically,  it  is  noted  that  the  spectrum  is  often  expressed  in  terms 
of  the  natural  frequency  v rather  than  the  angular  frequency  id.  Then, 


T (V) 
xx 


rxx  (,(v»  |*| 


r (v)  - 2tt  r (2ttv)  . 

xx  xx 


(5.24) 
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Altho  smoothing  is  required  to  reduce  the  variance  of  the  estimated  auto- 


spectrum, it  does  create  a distortion  of  the  estimate  ^^(co)  relative  to 


rxx(u>).  It  is  the  purpose  of  the  developments  of  this  section  to  quantify 


this  distortion. 


Two  weight  functions  are  frequently  used  in  smoothing:  one  due  to 


Bartlett — 


w (t)  = 1 - t/T  , 0 < t < T 
d m — — m 


t > T . 
in 


where  Tm  is  a lag  parameter,  and  one  due  to  Tukey — 


W (t)  = J (1  + COs(TTt/T  )),  0 < t <_  T 


t > T . 
m 


(5.26) 


In  the  following  analysis  the  simpler,  Bartlett  window  is  used.  However, 


the  effect  shown  on  the  estimated  autospectrum  is  representative  of  several 


lag  windows.  Also  in  this  analysis  we  take  the  estimated  autocovariance 


function  to  be  the  theoretical  autocovariance  of  a first-order  analog  process: 


Y (t)  = Yxx(0)  exp  (-Xt) . Specifically,  from  the  second  term  of  (2.11), 


the  autocovariance  of  the  scintillation  process  is 


Y (t)  - o2  e ~taBt 
XX  s 


(5.27) 


Then,  for  the  Bartlett  weight  function, the  smoothed  autospectrum  is 


, T 

2.  t m 


T (to) 

XX 


^ 1 

2qs 

* Jo 


cos(iot)  e Ust  (1  - t/x  ) dt. 

m 


(5.28) 


0 U T 

A 2 o 2 s ra 

r (w)  = i s 

XX  T-u£  Jo 


:os  (rr  x)  e"x  (1  - ) dx  . 


(5.29) 


From  (1.10),  the  constant  multiplier  on  the  r.h.s.  of  (5.29)  is  B/2tt. 
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I 


f fxx(“>  ■ f0  cos<^ x)  e'x  dx 


-J-f 

U)  T 

s m 1 0 


WT 
S IQ 


-X  \ . 

xe  cos( — x)  dx. 
00 

s 


After  some  manipulation. 


2tt 


where 


B rxB<“)  " rxx(“>  + “M 

— (jj  f 

+ e s m [a  (to)  cos  tot  + 8 (to)  sin  tot  1 , 

m m 


rxx(co)  = (i  + (<o/co8)2)-1 


-12  2 2 2-2 
a (to)  - toT  (to  - to)  (to  + to  ) 
sms  s 


a (to) 


1 - (to/to  )' 
s 


U)  T (1  + (to/to  )2)2 
s m s 


2-1  2 2-2 

8 (to)  = - 2 to  t to(to  + ao  ) 
s m s 


8(0))  = -2  (0)/0)  )(0J  T r1  (r  (to))2  . 

s s m xx 


(5.30) 


(5.31a) 


(5.31b) 


Note  that  the  normalized  smoothed  autospectrum  (2q/B)  T approaches 

the  unsmoothed  function  T as  T grows  infinite.  Some  numerical  examples 

xx  m 

A 

of  T were  evaluated  for  several  values  of  the  parameters  v (or  to  ) and 

s s 

T . Results  are  displayed  in  Figures  8 and  9.  The  effect  of  smoothing  is 
m 

to  shift  variance  (or  power)  to  lower  frequencies  and  to  introduce  a damped 

oscillatory  term.  These  distortions  are  not  very  significant  for  values 

of  v T greater  than  about  1.5. 
s m 
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Frequency,  V (hz) 

A 

Figure  8.  The  Smoothed  Autospectrum,  T (v)  , Due  to  a Bartlett  Lag  Window  Applied 
to  a First-Order  Autocovariance  With  Associated  Theoretical  (Unsmoothed)  Spectrum 
r<v),  (Parameters : v “ 1 hz,  t = 1 sec) 
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i 

i . i . 
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i . . , 
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. 

j 

Parameters 

Second-Ordt r Component 

The  digital  transfer  function  used  in  describing 
tracking  component  of  the  laser  spot  motion  (equation 
here. 


Hz(z) 


a 

o 


+ a. 


1 + b. 


-1 


+ b„ 


-2 

z 


the  second-order, 
(A. 5))  is  repeated 


(5 


The  autospectrum  of  the  process  is  proportional  to  the  squared  modulus 
|Hz(exp(jciiT))  |2.  Notationally , 

H2<»>  - |Hs(e>T)|2  - |H2(.J2”VI)|  2, 


with  sampling  generation  interval  T. 

Using  De  Moivre's  theorem, 

eJuT  m cos  u)T  + j s:n  ooT,  (5 

with  (5.32)  and  (5.33), 

2 2 
H.(ta)  - [(a  + a,  cos  U)T  + a.  cos  2ojT) 

2 O 1 2 

2 

+ (a^  sin  oil  + a^  sin  2toT)  ]/ 

( (1  + b^  cos  o)T  + b^  cos  2cjT)2  + (b^  sin  toT  + sin  2wT)2].  (5 


Using  equation  (4.6)  and  after  some  manipulation,  (5.35)  becomes 

4 

- w (3  + 4 cos  toT  + cos  2(jl)T ) 

H2(U)  = -iL_ ’ (5 

3 (to  + 1)  + 4(<jJ  - 1)  cos  coT  + (cj  + 1)  cos  2coT 

a a a 

where 

to  - tan  (irv  T) . (5 

cl  t 

Either  (5.35)  or  (5.36)  can  be  used  to  evaluate  the  squared  modulus  of 
the  second-order  dynamic  component.  Altho  (5.36)  has  a simpler  form 
than  (5.35),  computational  experience  using  single-precision  arithmetic 
on  the  IBM  360  series  computers  has  demonstrated  that  (5.36)  is  somewhat 
more  sensitive  to  truncation  error  than  (5.35).  In  fact,  it  is  noted 

(5-21) 
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.32) 


.33) 


.34) 


.35) 


.36) 


.37) 


that  generally  all  second-order  transfer  functions  and  squared  moduli! 
are  more  sensitive  to  loss  of  precision  due  to  arithmetic  truncation 
than  are  their  first-order  counterparts.  When  using  IBM  360  computers 
it  is  recommended  that  double-precision  arithmetic  be  used. 

The  above  expressions  for  the  second-orHer  squared  modulus  are 
normalized  so  that  h2(0)  = 1.  To  calculate  the  autospectrum,  ^xx(v)» 
of  the  second-order  (tracking)  process,  one  employs: 


rxx(v)  « H2(2itv)  Tnn(v)  , (5.38) 

where  r (v)  is  the  input  noise  spectrum  to  the  second-order  digital 

filter.  The  noise  spectrum  is  a function  of  the  sampling  generation  in- 

2 

terval,  T,  and  the  variance  of  the  tracking  process,  o^: 


T = 2T  a [y  (0)/y  (0)]  , 

nn  t nn  xx 


with  the  ratio  of  variances  y (0)/y  (0) 

nn  xx 


(5.39) 


given  by  (4.15a). 

As  with  the  first-order  process,  it  is  useful  to  make  numerical 
comparisons  between  the  squared  modulus  of  the  digital  implementation 
of  the  second-order  (tracking)  component  and  that  of  the  corresponding 
analog  process.  The  transfer  function  for  a general  second-order  analog 
process,  H (s) , is  given  by  equation  (4.1).  The  general  squared  modulus 
is,  then, 

H2(jco)  = 1/ [ (1  - (to/to  )2)2  + 4;2  (to/oo  )2]  (5.40) 

t a a 

With  the  j.itterworth  assumption  of  £ = 1//2  and  with  to  = 2ttv  , 

3 t 


H2(v)  = 1/(1  + (v/vt)4]  , (5.41) 

as  in  equation  (1.1)  with  A = 1. 

The  analog  gain  constant  A (equation  (1.5))  represents  the  analog 

input  noise  spectrum  and  is  comparable  with  r given  by  (5.39).  In  fact, 

lim  f “ A . (5.42) 

T* . r.  nn 


( 5-22) 


60 


Pursuing  the  example  of  Annex  1,  with 


I 


I 


\>t  *=  0.7  hz 
Ot  •=  0.23  ra 
T =0.1  sec, 

(j  ■ 0.219912,  and  from  (5.36): 

a 

(v)  = 2.33879  10~3 [3  + 4 cos(0. 62832V) 

+ cos(l. 25664v) ]/ [3.00716 
-3.99064  cos (0. 62832V)  + 1.002339  cos (1 . 25664v) ] . 


(5.43) 


Numerical  results  from  this  expression  and  similar  digital  squared  moduli! 
calculated  with  T = 0.05  sec  and  0.025  sec  are  compared  with  the  analog 
squared  modulus  in  Table  4.  Spectral  amplitudes  for  the  noise  are  also 
shown  in  Table  4.  Spectra  are  plotted  in  Figure  10. 


TABLE  4.  COMPARISON  OF  THE  ANALOG  AND  DIGITAL  SQUARED 
MODULI I OF  THE  SECOND-ORDER  (TRACKING)  ERROR  PROCESS 


V 

(hz) 

analog 

digital,  I 

^(v)  with 

T (sec) : 

H^(V) 

A=1 

0.025 

0.05 

0.10 

0.1 

1.000 

1.000 

1.000 

1.000 

0.2 

0.993 

0.993 

0.993 

0.994 

0.3 

0.967 

0.967 

0.968 

0.969 

0.4 

0.904 

0.904 

0.905 

0.907 

0.5 

0.793 

0.794 

0.794 

0.799 

0.6 

0.649 

0.650 

0.650 

0.653 

0.7 

0.500 

0.500 

0.500 

0.500 

0.8 

0.370 

0.369 

0.368 

0.365 

1.0 

0.194 

0.193 

0.191 

0.183 

1.2 

0.104 

0.103 

0.101 

0.092 

1.5 

0.0453 

0.0447 

0.0428 

0.0357 

2.0 

0.0148 

0.0144 

0.0132 

0.0089 

3.0 

0.00296 

0.00275 

0.00220 

0.00070 

5.0 

0.00038 

0.00031 

0.00015 

0.00000 

r .(v)  (m2/hz) 
nri 

0.0680 

0.0682 

0.0685 

0.0695 
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Autospectrum  (m  /hz) 


digital 


Frequency,  v (hz) 

Figure  10.  Autospectra  of  the  Second-Order  (Tracking)  Process(Parameters 
V = 0.7  hz,  0 = 0.230  m,  T = 0.1  sec) 
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MEMORANDUM  FOR  RECORD 


SUBJECT:  Cost  Differenli.il  (COSDIF)  Probnbilit ies 


1.  The  Systems  Analysis  1 ircctorate  was  requested  to  develop  ARMCOM' s 
Cost  Differential  (COSDtr)  Probabilities  for  both  provisioning  and  replen- 
ishment type  items.  This  MFR  addresses  the  COSDTF  probabilities  as  applied 
to  replenishment  icems. 

2.  The  COSDIF  formula  is  used  to  determine  if  it  is  cost  effective  to 
stock  an  item  when  the  expected  yearly  demand  is  12  or  less.  Inputs  to 
this  formula  include  various  costs  and  the  probability  that  an  item  will 
have  no  demands  in  the  next  two  years.  ARMCOM  is  presently  using  the 
DODI  4140.42  probabilities  (hard-coded  in  COSDIF)  which  were  developed 

by  the  DoD  from  data  supplied  from  one  Army  ICR  and  two  DSA  ICPs.  Beginnin 
in  January  1977,  ARMCOM  will  be  able  to  enter  probabilities  other  than 
those  identified  by  the  DODI. 

3.  In  1975,  TAR COM  recommended  a change  to  the  DoD  probabilities  based  on 
a study  of  TARCOM's  managed  items,  TARCOM  Report  No,  75-35,  subject: 

COSDIF  Probabilities  for  TARCOM.  Table  1 presents  th.j  probabilities  as 
identified  in  DODI  4140.42  and  those  proposed  by  TARCOM.  The  proposed 
probabilities  were  developed  from  a study  of  four  years  of  demand  history 
records  on  approximately  2,000  TARCOM  managed  items.  TARCOM's  probabilitie 
are  considerably  lower  than  those  arrived  at  by  the  DoD  and  are  now  being 
evaluated  by  TARCOM  as  to  their  affect  on  the  number  of  items  stocked, 
budget,  readiness  ratio  and  etc, 

4.  The  approach  used  in  developing  ARMCOM* s probabilities  for  replenish- 
ment items  was  to  obtain  the  demand  history  of  ARMCOM  managed  items. 
However,  the  only  accessible  data  was  the  Nov  76  ALPHA  Demand  Return  Dis- 
posal (DRD)  file,  which  contains  only  2 years  of  demand  data  (from  Dec  74 
to  Dec  76)  for  the  9,000  ARMCOM  managed  active  items.  Since  the  ALPHA 
system  keeps  only  the  current  2 years  of  data,  a search  was  made,  within 
the  Command,  to  obtain  the  oldest  Operations  Readiness  Oriented  Supply 

St  itus  (OROSS)  reports.  A Jan  76  report  was  located  which  summarizes 
demands  over  a two  year  pei rod  (Jan  74  through  Jan  76)  and  if  merged 
into  the  DRD  files,  will  provide  a 3 year  data  base.  Since  only  the  DRD 
files  are  accessible  at  this  tine,  the  probabilities  as  developed  in  this 
MFR  are  based  on  only  the  two  y( ars  of  demand  data. 

5.  The  assumptions  and  criteria  used  in  developing  the  replenishment 
probabilitie:-,  are  as  follows: 
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a.  The  average  yearly  demand  (demand  frequency)  of  an  item  is 

estimated  to  be  the  number  of  demands  (requisitions)  that  occurred  in  the  first  of 
the  two  years  of  data,  i.e,,  Dec  74  through  Dec  75. 

b.  The  demand  probabilities  remain  constant  for  the  application 
period,  approximately  2 years. 

6.  The  ARMCOM  probabilities  were  developed  by  determining  for  each  group 
of  item:,  having  the  same  demand  frequency,  the  proportion  of  the  group 
having  no  demands  in  the  succeeding  year.  For  example,  if  100  items  had 
a demand  frequency  of  6 in  the  first  year  and  40  of  those  items  had  no 
demand  in  the  2nd  year,  the  proportion  would  be  0.4.  This  proportion  is 
an  estimate  of  the  probability  of  having  no  demands  for  one  year  given 
the  demand  frequency  for  the  prior  year.  The  square  of  this  proportion  is 
then  an  estimate  of  the  probability  of  no  demands  in  two  years.  Table  2 
presents  a summary  of  ARMCOM' s 2 years  of  demand  history;  the  number  of 
those  items  having  a demand  frequency  of  1 through  12  in  the  first  year; 

the  number  of  those  items  having  no  demand  the  2nd  year;  and,  the  probability 
of  no  demands  in  the  next  succeeding  years.  As  shown  in  Table  2,  an  item 
ha'ing  a demand  frequency  of  1 has  a probability  0.58  of  having  no  demands 
over  a one  year  peri  id  and  a probability  of  0.34  of  having  no  demands  over 
a two  year  period.  As  den. and  frequency  increases  the  probability  of  no 
demands  decreases  sharply. 

7.  There  are  three  options  available  to  ARMCOM  (lncl  4): 

a.  Continue  to  use  the  D0DI  4140.42  probabilities. 

b.  Implement  the  probabilities  developed  by  TARC0M. 

c.  Implement  the  proposed  ARMCOM  probabilities. 

Table  3 presents  each  option.  There  is  a significant  difference  between 
the  D0DI  set  and  those  listed  under  1ARC0M  and  ARMCOM,  the  latter  two  being 
considerably  lower  than  DoD's,  However,  they  do  not  appear  to  be  totally 
unreasonable,  DA  Pamphlet  70-5,  'Mathematics  of  Military  Action,  Operations 
and  Systems",  states,  "...  an  item  whose  average  demand  rate  is  one  piece 
per  time  period  will  have  a probability  of  nearly  0.4  of  showing  a history 
of  no  d cm.  nd  at  all  during  another  time  period  of  the  same  length,"  Although 
the  stater ent  refers  to  quantities,  an  item  that  has  a demand  frequency  of  1 per 
period  will  have  at  least  an  average  demand  rate  of  2 piece  per  period.  A higher 
piece  demand  rate  would  logically  lower  the  probability  of  no  demand; 

ARMCOM'.';  average  number  of  pieces  per  demand  incidence  is  38.  The  probability 
of  no  demand  in  a two  year  period  for  an  item  having  a demand  frequency  of  1 
is  0.44  for  TAPOPM  item:;  and  0.3';  for  ARMCOM  items  which  compares  favorably 
with  the  probability  of  0.  i as  staled  in  DA  Pamphlet  70-5.  In  comparing  the 
proposed  ARMCOM  and  fAROOM  probabilities;  ARMCOM 's  are  lower  and  converge 
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to  zero  very  sharply  (having  a value  of  zero  for  demand  frequencies  greater 
than  7)  where  TARCOM's  probabilities  remain  relatively  constant  ranging  from 
0.04  to  0.11  for  demand  frequencies  of  4 to  12.  These  differences  are  caused 
by  various  factors,  e.g.,  items  peculiar  to  each  Command,  degree  of  item 
commonality,  or  total  number  of  items  stocked.  However,  the  impact  on  ARMCOM 
in  either  implementing  the  proposed  TARCOM  or  the  proposed  ARMCOM  probabili- 
ties would  certainl'-  be  an  increase  in  the  number  of  items  stocked  which  in 
turn,  affects  ARMCOM's  budget,  and  readiness  posture.  The  degree  of  the 
impact  and  the  extent  of  the  impact  in  dollar  savings,  to  the  Army  however, 
is  dependent  on  the  accuracy  of  the  probabilities  in  predicting  no  demands. 

8.  In  summary,  the  probabilities  prepared  by  TARCOM  and  those  developed 
for  ARMCOM  agree  in  that  b 5th  are  considerably  lower  than  those  identified 

in  the  D0D1  4140.42.  However,  ARMCOM's  probabilities  do  differ  from  TARCOM's; 
they  are  generally  lower  aid  reach  a value  of  0.0  for  a demand  frequency 
of  8 or  more  (see  Table  3),  but  these  differences  may  be  representative  of 
the  differences  for  each  Cimmand's  managed  items. 

9.  Recommendations  and  Conclusions: 

a.  A study  be  conducted  to  determine  the  dollar  impact  of  implemencing 
the  proposed  ARMCOM  probabilities. 


b.  The  proposed  ARMCOM  probability  set  be  updated  every  year  for 
approximately  3 to  5 years  and  then  reviewed  every  2 years,  thereafter. 
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TABLE  1.  COSPIF  PROBABILITIES 


DEMAND 

FREQUENCY*  1)01)1  41 4 0.4 2 PROPOSED  BY  TARCOM 


1 

.82 

.44 

2 

.67 

.32 

3 

.55 

.21 

4 

.45 

.11 

5 

.37 

.10 

6 

.30 

.10 

7 

.25 

.07 

8 

.20 

.08 

9 

.17 

.04 

10 

.14 

.05 

11 

.11 

.04 

12 

.09 

.06 

*Only  those  items  with  demnnd  frequency  of  12  or  less  were 
considered . 
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TABU:  3. 

OORniF  PRO  BA 

;n.mus  for  repleeio;/::: 

N'T  type  it;;-:;; 

1 

1 

DEMAND 

l'REOBEXCV-- 

ROOT.  A 14 0,42 

PROPOSED  BY 
TARCOM 

pro; '03!  j 
A:\y.ar.  - v 

* 

i 

1 

.32 

.44 

.34 

t 

2 

.67 

.32 

.15 

\ 

3 

.35 

.21 

.05 

4 

.45 

.11 

.04 

5 

.37 

.10 

.02 

\ 

6 

.30 

.10 

.01 

7 

.25 

.07 

.01 

8 

.20 

.08 

0,0 

9 

.17 

.04 

0.0 

10 

.14 

.05 

0.0 

11 

.11 

.04 

0.0 

12 

.09 

.06 

0.0 

-•‘Only  those  items 

with  demand 

frequency  of  12  or  less 

were  considered. 
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SUBJECT:  Change  to  Cost  Differential  (COSDIF)  Probabilities 


SHE  DISTRIBUTION 


| 1 . References : 

a.  TARCOM  Report  No.  75-35,  subject:  COSDIF  Probabilities  for 

TARCOM,  attached  as  Inclosure  1. 

b.  System  Change  Request  #XSUZZC528701 , subject:  COSDIF 

Probabilities,  14  October  1975,  attached  as  Inclosure  2. 

c.  AMCSU-KP  letter,  5 December  1975,  subject:  SCR  for  Conditional 

Probabilities  for  Provisioning  Items,  attached  as  Inclosure  3. 

2.  The  purpose  of  this  letter  is  to  provide  information  regarding  a 
potentially  significant  change  to  the  COSDIF  stockage  model  and  provide 
guidance  on  the  use  oi  this  change.  At  present,  one  s>.t  of  COSDIF 
probabilities  is  applied  to  both  provisioning  items  and  replenishment 
type  items.  This  set  of  probabilities,  which  is  hard-coded  in  COSDIF, 
is  identified  in  DODI  4140.42. 

3.  Beginning  in  January  1977,  CCSS  will  be  able  to  accommodate  two 
sets  of  COSDIF  probabilities.  One  set  will  apply  to  provisioning  items 
and  will  remain  hard-coded  until  April  1977  (Release  49).  The  other 
set  will  apply  to  replenishment  type  items  and  can  be  entered  by  the 
NICP  following  release  47  in  January  1977. 

4.  In  applying  this  added  capability,  there  are  several  options  open 
to  the  commodity  commands.  The  NICP  may: 

; a.  Continue  to  use  the  set  of  probabilities  cited  in  DODI  4140.42 

for  both  categories. 

b.  Develop  their  own  probabilities, 

c.  Elect  to  use  for  provisioning  items  the  set  developed  for  the 
Navy  by  ALRAND.  In  no  case  should  this  set  be  used  for  rcplcnishn cr.t 
i terns . 
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d.  Elect  to  use  the  set  developed  by  TARCOM  for  replenishment 
items.  In  no  case  should  this  set  be  used  for  provisioning  items. 

5.  In  all  cases  where  the  command  elects  to  use  probabilities  other 
than  those  identified  in  the  DODI,  the  command  should  make  an  analysis 
to  insure  probabilities  selected  are  appropriate  and  assess  budget 
impact.  The  Inventory  Research  Office  and  ALMSA  will  provide  assistance 
in  this  analysis.  The  command  should  also  be  prepared  to  provide 
justification  for  such  selections,  should  the  question  arise  during 
budget  hearings. 

6.  Any  questions  regarding  this  guidance  should  be  discussed  with 

LTC  D.  Weir,  DRCMM-RS,  AUTOVON  284-8221.  ; 

FOR  THE  COMMANDER: 


3 Incl  ' ' 
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v/c d r^U S_ .A rmyj\ rma men.t s ...Comma qcU„ATT >{ RSAf^JiM-P ,„Ro c.LJUla nd.  IL  61201 
Cdr,  US  Army  Aviation  Systems  Command,  ATTN:  DRSAV-QP,  St.  Louis, 

MO  63166 

Cdr,  US  Army  Electronics  Command,  ATTN:  DRSEL-MM-P,  Ft  Monmouth,  NJ  07703 

Cdr,  US  Army  Missile  Command,  ATTN:  DRSMI-SOT,  Redstone  Arsenal, 

AL  35809 

Cdr,  US  Army  Tan <-Automotive  Readiness  Command,  Warren,  MI  40090,  ATTN:  DRSTA- 
Cdr,  US  Army  Troop  Support  Conmand,  ATTN:  DRSTA-SX,  4300  Goodfellow  Blvd, 

St.  Louis,  MO  63101 

Cdr,  US  DARCOM  Automated  Logistics  Systems  Agency,  ATTN:  DRXAL-MCC,  St. 

Louis,  MO  63101 

Ch,  DRC  Inventory  Research  Office,  Philadelphia,  PA  19106 


Colonel , GS 
Assistant  for 
Secondary  Items  Programs 
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Commander 

US  Army  Materiel  Development  and  Readiness  Command 
ATTN:  DRCMA 

DRCRD 
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Commander 
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ATTN:  DRSAR-CG 
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DRSAR-EN 
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DRSAR-PA 
DRSAR-PP 
DRSAR-OP 
DRSAR-QA 
DRSAR-IS 
DRSAR-MM 
DRSAR-MA 
DRSAR-AS 
DRSAR-SF 
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Commander 
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US  Army  Electronics  Command 
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Commander 
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P.  0.  Box  209,  Main  Office 

St.  Louis,  MO  64502 

1 Commander 

US  Army  Troop  Support  Command 

ATTN:  DRSTS-G 

St.  Louis,  M0  63120 

1 Project  Manager  fo  - Cannon  Artillery  Weapons  Systems 

ATTN:  DRCPM-CAWS 

Dover,  NJ  07801 

1 Commander 

US  Army  Development  and  Readiness  Command 

Office  of  the  Project  Manager  for  Selected  Ammunition 

ATTN:  DRCPM-SA 
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ATTN : DRCPM-PBM 
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1 Product  Manager  for  Advanced  Attack  Helicopter  Systems 
US  Army  Aviation  Systems  Command 
St.  Louis,  M0  63166 

1 Product  Manager  for  AH-I  Cobra  Series  Aircraft 
US  Army  Development  & Readiness  Command 
P.  O.  Box  209 
St.  Louis,  M0  63166 
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1 Commandant 
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1 Commandant 
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1 Commander 

US  Army  Missile  & Munitions  Center  & School 
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1 Director 
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1 Commander 

US  Army  TRADOC  Systems  Analysis  Activity 
White  Sands  Missile  Range 
White  Sands,  NM  88002 

1 Director 

Advanced  Research  Projects  Agency 
1400  Wils on  Boulevard 
Arlington,  VA  22209 

1 Commander 

Defense  Logistics  Studies  Information  Exchange 
Fort  Lee,  VA  238C  L 

1 Commander 

US  Army  Logistics  Management  Center 
ATTN:  DRXMC-LS 

Fort  Lee,  VA  23801 

1 Commander 

US  Army  Logistics  Center 

ATTN:  ATCL-S 

Fort  Lee,  VA  2380L 

12  Defense  Documentation  Center 

Cameron  Station 
Alexandr  .a,  VA  22  314 
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